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Summary. We study the motion of electrons in a periodic background potential 
(usually resulting from a crystalline solid). For small velocities one would use either 
the non-magnetic or the magnetic Bloch hamiltonian, while in the relativistic regime 
one would use the Dirac equation with a periodic potential. The dynamics, with the 
background potential included, is perturbed either through slowly varying external 
electromagnetic potentials or through a slow deformation of the crystal. In either 
case we discuss how the Hilbert space of states decouples into almost invariant 
subspaces and explain the effective dynamics within such a subspace. 



1 Introduction 

In a crystalline solid the conduction electrons move in the potential created 
by the ions and the core electrons. Somewhat mysteriously and linked to the 
Pauli exclusion principle, the Coulomb repulsion between conduction electrons 
may be ignored, within a good approximation. Thereby one arrives at a fun- 
damental model of solid state physics, namely an ideal Fermi gas of electrons 
subject to a periodic crystal potential. Let r be the periodicity lattice. It is 
a Bravais lattice and generated through the basis {71,72,73}, 7j € M 3 , as 

3 

r = {7 = a ni with aez 3 }. (1.1) 

The crystal potential Vr is then /^-periodic, i.e., Vr ■ R 3 — ► K and Vr(x + 
7) = Vr{x) for all 7 6 r, and the electrons are governed by the one-particle 
hamiltonian 

H SB = ~A X + V r . (1.2) 

Hqb is the (Schrodinger)-Bloch hamiltonian. A wave function ip t G L 2 (M. 3 ) 
evolves in time according to the Schrodinger equation 
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We have chosen units such that the mass of an electron m c = 1 and h = 1. 
The electron charge, e, is absorbed in Vr- Since Vr is periodic, electrons move 
ballistically with an effective dispersion relation given by the Bloch energy 
bands E n , see below for a precise definition. E n is periodic with respect to 
the lattice r* dual to r, E n (k + 7*) = E n (k) for all 7* e T*, k e M 3 . This 
feature makes the dynamical properties of a Bloch electron very different from 
a massive particle with dispersion Ef iee (k) = \k 2 valid in case Vr = 0. 

The thermodynamics of the electron gas is studied taking iJge as a starting 
point. Dynamically, however, one wants to probe the response of the electrons 
to external forces which very crudely come in two varieties. 

(i) External electromagnetic potentials. Electrostatic potentials manufactured 
in a lab have a slow variation on the scale of the lattice r. Therefore we 
set T4xt(a0 = e<j)(ex), e the charge of the electron, with e a dimensionless 
parameter and cj> independent of e. e <C 1 means that the potential V ex t has a 
slow variation when measured with respect to the lattice spacing of r. Note 
that the electrostatic force is 0(e) and thus weak. External magnetic fields 
on the other hand can be so strong that the radius of gyration is comparable 
to the lattice spacing. It then makes sense to split the vector potential as 
A n + A CKt , where A (x) = — ^B A x with B S M 3 a constant magnetic field. 
Included in -ffge, this yields the magnetic Bloch hamiltonian 

Hub = A ) 2 + V r . (1.4) 

A cxt is a probing vector potential in addition to Aq. A cxt is slowly varying 
on the scale of the lattice, A ext (x) = A(ex) with A independent of e, and the 
corresponding magnetic field is small of order e. Including all electromagnetic 
potentials, for simplicity with the electric charge absorbed into A and <fi, the 
hamiltonian becomes 

H = ^(-iV x - A {x)- A(ex)) 2 + V r (x)+<j)(ex). (1.5) 

(ii) Mechanical forces. The crystal lattice can be deformed through exter- 
nal pressure and shear. Thereby an electric polarization is induced, an ef- 
fect which is known as piezoelectricity. If charges are allowed to flow, in this 
way mechanical pressure can be transformed into electric currents. The me- 
chanical forces are time-dependent but slow on the typical time-scale of the 
electrons. Therefore in (|1.2p Vr(x) is replaced by Vr( e t) (x, et). r(et) is the 
instantaneous periodicity lattice and is defined as in (jTTTj) . Vp(t) is space- 
periodic, i.e. Vr(t)(x + 7,£) = V(x,t) for all 7 G r(t). The special case of a 
time-independent lattice, r(t) = _T, but a still slowly in time varying crystal 
potential is also of interest. For example, one may imagine a unit cell with 



Periodic structures 3 



two nuclei. If the two nuclei are displaced relative to each other, then r re- 
mains fixed while the crystal potential in the unit cell changes with time. The 
resulting piezoelectric hamiltonian reads 

H PE (t) = --A x + V r(st) (x > et). (1.6) 

Our general goal is to understand, in each case, the structure of the solu- 
tion of the time-dependent Schrodinger equation for small e. Obviously, H in 
(|1.5|> is a space-adiabatic problem, while (II. 6p corresponds to a time-adiabatic 
problem. However in the latter case it turns out to be profitable to transform 
to a time- independent lattice, say -T(O). Then also terms varying slowly in 
space are generated. Thus, in the general case the full power of the space- 
adiabatic perturbation theory PST03a, ITeu03| will be needed. A word of 
caution must be issued here for the magnetic Bloch hamiltonian. To use the 
methods from [PST03a] in this context, the unperturbed Hamiltonian must 
be periodic, which is the case only if the magnetic flux per unit cell is rational. 
One can then define an enlarged magnetic unit cell such that Hmb is invariant 
with respect to magnetic translations. If the magnetic flux is not rational, the 
crutch is to include in Aq a nearby rational flux part of the magnetic field, 
with a small denominator, and to treat the remainder as A ext - 

To achieve our goal, depending on the context we use one of the peri- 
odic hamiltonians as backbone. The periodic hamiltonian is denoted by H per 
with ffper either Hsb, or Hmb, or TJpe at fixed t, or TJls from (|1.8[) . or iJrjB 
from (|1.9p . As explained below, the Hilbert space TC = L 2 (R 3 ) then splits 
as Ti. = (B^Lo^™' where- n is the band index. Each subspace 7i n is invari- 
ant under exp[— itH per ] and H pcr restricted to TL n is unitarily equivalent to 
multiplication by E n (k). E n (k) is the effective hamiltonian associated to the 
u-th band. The complexity of the full problem has been reduced substan- 
tially, since only a single band dynamics has still to be studied. Modifying 
flper such that it becomes slowly varying in space-time is, vaguely speaking, 
a small perturbation. Thus one would expect that the invariant subspace 7i n 
is to be substituted by a slightly tilted subspace. On this subspace E n (k) 
will turn into a more complicated effective hamiltonian. The difficulty is that, 
while the dynamics generated by H pel can be computed by solving a purely 
spectral problem, none of the perturbed hamiltonians can be understood in 
this way. In particular, one has to spell out carefully over which time scale 
the slightly tilted subspace associated to TC n remains approximately invari- 
ant and in what sense the dynamics generated by the effective hamiltonian 
approximates the true time evolution. 

To lowest order the effective hamiltonian can be guessed from elementary 
considerations and belongs to a standard tool of solid state physics |AM76| . 
The guess provides however little hint on the validity of the approximation. 
There one needs a mathematical theorem which states precise conditions on 
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the initial wave function and provides an error bound, from which the time 
scale for validity can be read off. 

Under the header "geometric phase" physicists and quantum chemists have 
realized over the past twenty years, say, that the first order correction to the 
effective hamiltonian carries a lot of interesting physics, see [BMKNZ03 for a 
recent comprehensive overview. For the magnetic Bloch hamiltonian the first 
order correction yields a Hall current proportional to the Chern number of the 
magnetic Bloch vector bundle. Similarly the modern theory of piezoelectricity, 
expresses the piezocurrent as an integral of the Berry connection over the 
Brillouin zone, see King-Smith, Vanderbilt [KSV93] and Resta [Res94j . First 
order effective Hamiltonians are no longer guessed so easily and it is convenient 
to have the systematic scheme }PST03aj available. 

In nature electrons are spin ^ particles. The wave function is thus (de- 
valued and the hamiltonian in (| 1 . 5|) is modified to 

H = i ( - iV x - A (x) - A{ex)f + V r {x) + <f>(ex) -^a-(B Q + eB{ex)) (1.7) 

with B — V A A for the slowly varying part of the magnetic field. Here 
a = ((7i, (i2, 03) is the 3-vector of Pauli spin matrices. Besides the term pro- 
portional to the uniform magnetic field Bo, H acquires a subleading term of 
order e. More accurately one may want to include the spin-orbit coupling. 
The periodic piece of the hamiltonian reads then 

#ls = --A* + V r (x) + ~a • (VVr(z) A (-iV*)) (1.8) 

and the slowly varying potential is added as in (jTTTJ) with the additional sub- 
leading term e\a ■ (V(p(ex) A (—iV x ))- 

Depending on the crystalline solid, the conduction electrons can move so 
fast that relativistic corrections become important. On the one-particle level 
an obvious choice is then the Dirac equation with a periodic potential Vr- 
Wave functions are C 4 -valued and the hamiltonian reads 

#db = I3m c c 2 + ca ■ p + V r , p=—iV x . (1.9) 

We introduced here the mass, m , of the electron and the speed of light, c. 
The 4x4 matrices /3, ai, 0:2, 0:3 are standard and defined in |Tha94( lYnd96] . 
for example. Note that the Lorentz frame is fixed by the solid, i.e. by Vr- 

In fact, the non-relativistic limit for TJdb yields the spin-orbit hamiltonian 
Hls |Tha94[ lYnd96j . If ||Vr|| is bounded, for sufficiently large c, the Dirac 
hamiltonian i?DB has a spectral gap, which widens as c — > 00. Projecting 
onto the electron subspace, to leading order in 1/c one obtains the Pauli- 
Bloch hamiltonian — (l/2m e )A x + Vr with the spin-orbit coupling in (|1.8j) as 
a correction of strength 1/(to c c) 2 . In addition the crystal potential is corrected 
by -A x Vr(x)/8(m c c) 2 . 

In our contribution we will provide some background on how to establish, 
including error bounds, the validity of the approximate dynamics as generated 
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by an effective hamiltonian, including order e corrections, for most of the 
models mentioned in the introduction. For this purpose it is necessary to 
briefly recall the spectral theory for the periodic hamiltonian, which is done 
in the following section. In the subsequent sections we deal with particular 
cases in more detail. We start with the non-magnetic Bloch hamiltonian, see 
(|1.5jl with Bq = 0. For the magnetic Bloch hamiltonian we explain how Bq — > 
and Bo — > oo may be viewed as particular adiabatic limits. Piezoelectricity is 
discussed in the last section. 

Remark. Our contribution is one part of the research project jointly with 
S. Bauer and M. Kunze within the Schwerpunkt. Their part will be covered in 
[BK06]. Both contributions appear now as almost disjoint, which only reflects 
that we wanted to present a coherent story. The unifying aspect is an adiabatic 
limit for wave-type evolution equations. In this contribution we stay on the 
level of effective hamiltonians while in [BK06 one pushes the scheme to the 
first dissipative correction. 

2 The periodic hamiltonians 

We consider a general dimension, d, and assume that the periodicity lattice 
r is represented as 

r = jx e M. d : x = Ej=i«i 7j for some a G Z d | , (2.1) 

where {71, . . . , 7^} are vectors spanning K d . We denote by r* the dual lattice 
of r with respect to the standard inner product in R d , i.e. the lattice generated 
by the dual basis ... ,7^} determined through the conditions 7* • jj — 
2n5ij, i,j € {1, . . . ,d}. The centered fundamental domain M of r is defined 

by " ' 

M = {x€R d :x = Y£=\Pi 1j ^ ft G [~h 3]} > ( 2 - 2 ) 

and analogously the centered fundamental domain M* of r* . The set M* is 
the first Brillouin zone in the physics parlance. 

Assumption 1. The crystal potential Vr '■ R — > K satisfies Vr(x + 7) = 
Vr(x) for all j G r , x £ M. d . Vr is infinitesimally bounded with respect to 
-A. 

It follows from Assumption 1 that the periodic hamiltonians discussed 
below are self-adjoint on the domain of —A. 

2.1 The Bloch hamiltonian 



We consider 



H = --A + V r . 



(2.3) 
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The periodicity of H is exploited through the Bloch-Floquet-Zak transform, 
or just the Zak transform for sake of brevity |Zak68j . The advantage of such 
a variant is that the fiber at k of the transformed Hamiltonian operator has 
a domain which does not depend on k. 
The Zak transform is defined as 

(WzVO(M) := 5>- ifc -(*+^0c + 7), (k,x)eR 2d , (2.4) 

initially for a fast-decreasing function ip £ S(M. d ). One directly reads off from 
(|2.4[) the following periodicity properties 

(Uzif>)(k,y + i) = (Uzil>)(k,v) for all 7 ef, (2.5) 

(U z tp)(k + X,y) =e- iy - x (U z ^)(k,y) for all X E T* . (2.6) 

From (|2.5p it follows that, for any fixed k S R d , \lAz'ip)(k,-) is a /'-periodic 
function and can thus be regarded as an element of Hf = L 2 (M). M — 
M. d /r and it has the topology of the d-dimensional torus T d . On the other 
side, Equation (|2.6| involves a unitary representation of the group of lattice 
translations on _T* (isomorphic to r* and denoted as A), given by 

r-.A^UCHt), A^r(A), (r(X)<p) (y) = e^Mv)- (2-7) 
It is then convenient to introduce the Hilbert space 

H T = {i> S L 2 oc {R d , Hi) : ip(k - X) = r(A) tp{k) for all A e A} 

= L 2 T (R d , H f ) , (2.8) 

equipped with the inner product 

(i>, <p)h t = f dk(Tjj(k), <p(k)) Ht . (2.9) 

Obviously, there is a natural isomorphism between 7i T and L 2 (M*, Tt{) given 
by restriction from K d to M*, and with inverse given by r-equivariant con- 
tinuation, as suggested by (|2.6p . Equipped with these definitions, one checks 
that the map in (|2.4p extends to a unitary operator 

U z : L 2 (R d ) ^H T = L 2 (M*,L 2 (M)), (2.10) 

with inverse given by 

(WzV)(aO = / dke ik -^(k,[x}), (2.11) 
Jm* 

where [•] refers to the a.c. unique decomposition x = j x + [x], with ^ x G r 
and [x] € M. 
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As already mentioned, the advantage of this construction is that the trans- 
formed hamiltonian is a fibered operator over M* . Indeed, for the Zak trans- 
form of the hamiltonian operator (|2.3p one finds 

r® 

UzHUz 1 = \ dkH(k) (2.12) 

with fiber operator 

H(k) = ^(-iV y + k) 2 + V r (y), k e M* . (2.13) 

By Assumption 1, for fixed k S M*, the operator H(k) acts on L 2 (M) 
with the Sobolev space H 2 (M) as domain independently of k 6 M*. Each 
fiber operator H(k) has pure point spectrum accumulating at infinity. For 
dcfinitcncss the eigenvalues are enumerated according to their magnitude 
Eo(k) < Ei(k) < E2(k) < ... and repeated according to their multiplicity. 
E n : M* — > R is the n-th energy band function. It is continuous on M* when 
viewed as a d-torus. Generically the eigenvalues E n {k) are non-degenerate. 
Of course, there may be particular points in fc-space where particular energy 
bands touch each other and the corresponding eigenvalue becomes degenerate. 
The normalized eigenfunction corresponding to E n (k) is the Bloch function 
and denoted by tpn(k) £ H 2 (M). It is determined only up to a fc-dependent 
phase factor. A further arbitrariness comes from points where energy bands 
touch. We denote by P n {k) the projection along ip n (k) and set 

P n = dkP(k), H n = P n L 2 (R d ). (2.14) 

Through the Zak transform we have achieved the product structure 

H = H S ®H{, H s = L 2 {M*), H{ = L 2 (M). (2.15) 

ijj G Tin is of the form 4>(k)ip n (k, y). The band index n fixes the local pattern of 
the wave function ip while 4>(k) provides the slow variation. Therefore L 2 {M*) 
is the Hilbert space of the slow degrees of freedom. On the other hand for 
fixed k, one has oscillations in time determined by the eigenvalues E n (k). On 
long time scales, these become fast oscillations. Therefore Tif = L 2 (M) is the 
Hilbert space of the fast degrees of freedom. 

Since [P n ,H] = 0, the subspaces TL n are invariant under e~ lHt . P n e~ lHt P n 
is unitarily equivalent to multiplication by e -15 '™^) 4 on L 2 (M*). Note that, in 
general, 7i n is not a spectral subspace for H . The band functions generically 
have overlapping ranges. Therefore, if slowly varying terms are added to the 
hamiltonian, the dynamics can no longer be captured so easily by a spectral 
analysis of the perturbed hamiltonian. 
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2.2 The magnetic Bloch hamiltonian 

We consider d = 3. The hamiltonian reads 

H = ^(-iW x -A(x)) 2 + V r (x), xeR 3 , (2.16) 

with A(x) — —\B Ax, Bel 3 . Physically the most relevant case is d = 2. It 
is included here by setting x = (xi,x 2 ,0) and B = (0,0, B Q ). Following Zak 
|Zak64j . see also [DGR02] , one introduces the magnetic translations 

(T a il>)(x) = (e- ia <- iV - +A ^^)(x) = e laA{x H{x - a) (2.17) 

with a£l 3 . They satisfy the Weyl relations 

T a Tn = e-^ B < aA ^T a+l3 = e- iB < ah ^T p T a . (2.18) 

To have a commuting subfamily we need 

Assumption 2. The magnetic field B is such that B ■ (7 A 7') G 2ttQ for all 

7,7' er. 

In the two-dimensional case our assumption requires that the magnetic 
flux per unit cell, Bq • (71 A 72), is a rational multiple of 2ir. 

Under the Assumption 2 there exists a sublattice To C T such that B ■ (7 A 
7') G 2irZ for every 7, 7' G To- To is not unique. The set {T a } a er is a family 
of commuting operators, which commute with H. Since T a Tp = ±T a +/3, the 
magnetic translations still form only a projective group. It becomes a group 
by an even smaller sublattice A c To such that B ■ (7 A 7') G AirZ for all 
7,7' G Ti. Another common choice is a further modification of the phase 
through 

T a = e-if^Ta (2.19) 

with <p(a) = BYCt20Lz+B?,OL\a2 — B20i\(Xz - Then T a Tp — T a+ p for all a,/3 G To. 

We can now proceed as in the non-magnetic case. The Zak transform 
becomes 

(2/zV0(fc, x ) = z~ tk ' [x+l) T~fHx) , (k, x) G R 6 . (2.20) 

76-fo 

The properties of Uz^ are as in (|2.2p . (|2.3p provided T is replaced by T , 
and H T is replaced by H B = {u G if oc (R d ,H f ) : (|2~23|) below holds true}. In 
particular, H of (12.16|) admits the fiber decomposition 

UzHUz 1 = / dkH{k) (2.21) 
with M* the first Brillouin zone of Tq and with the fiber operator 

H(k) = ~(-%V y + \b A y + k) 2 + V r {y) ■ (2.22) 
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The domain of H(k) is independent of k but, in contrast to H(k) from (|2.13[) . 
a function u in the domain has to satisfy the more complicated boundary 
condition 

e-i y < aAB ^u(y-a)=u(y). (2.23) 

2.3 Dirac hamiltonian, spin-orbit coupling 

The Dirac hamiltonian with periodic potential reads 

H = j3-ia-V x + V r (x), (2.24) 

where we have set m c = 1, c = 1. As for the Bloch hamiltonian, H admits the 
fiber decompositon 



H= / dkH(k) (2.25) 



with fiber hamiltonian 



H(k)=/3 + a-{-iV y + k) + Vr(y). (2.26) 

H(k) acts on L 2 (M, C 4 ) with periodic boundary conditions (|2.5p . The free 
Dirac operator has a spectral gap of size 2, in our units, between the elec- 
tron and positron subspace. If we assume ||Vr|| < 1, then this gap persists 
and the eigenvalues can be labelled as Eo(k) < Ei(k) < . . . in the electron 
subspace and as E-\(k) > E-2(k) > ... in the positron subspace. One has 
E_i(k) < Eo(k) for all k £ M*. (In fact the labelling can be achieved without 
a restriction on ||Vr||, see [Mau03 ). 

For Vr — 0, the eigenvalue E(k) is two-fold degenerate reflecting the spin 
i of the electron, resp. positron. This degeneracy persists if the periodic po- 
tential is inversion symmetric, see IMau03j for details. 

Proposition 2.1. Let H be given by &2.24\ ) with \\Vr\\ < 1. Let there exist 
a G R 3 such that Vr{x + a) = Vr{— x + a). Then each E n (k) is at least 
two-fold degenerate. 

Proof. Without loss of generality we may assume a — 0. We use the standard 
basis for the a-matrices, see |Th a94 . In this basis time-reversal symmetry is 
implemented by the anti-unitary operator 

Tip(y) = -ia 3 a 1 ip*(y), (2.27) 

where the complex conjugation is understood component-wise. Using that 
aia^ai = —a^aiae, I = 1,2,3, where ai refers to matrix element-wise com- 
plex conjugation, one checks that 

— i\7 y aiT — —iT^yOLi , kaiT = —Tkag (2.28) 

and therefore 
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T- x H(k)T = H(-k) . (2.29) 
Secondly we use space inversion as 

Ri>{y) = m-v) ■ (2-30) 

Then 

R- l H{k)R = H(-k) . (2.31) 
Combining both symmetries implies 

T~ 1 R- 1 H(k)RT = H(k) . (2.32) 

If H(k)tp = Eip, then also RTip is an eigenfunction with the same eigen- 
value. Thus our claim follows from (ip, RTip) = 0. To verify this identity we 
note that — ia^ai = diag (<72)C2) and (\, R&2X*) = f° r an arbitrary two- 
spinor X- 

Corollary 2.2. The eigenvalue E n (0) of H(fi) is at least two-fold degenerate. 

Proof. Since T*H(0)T = H(0) by (|2~29f and (<tp, (-ia 3 ai)ip*) = 0, the claim 
follows. 

If Vr is not inversion symmetric, generically an energy band is two-fold de- 
generate at k — and then splits into two non-degenerate bands. Note that a 
non-degenerate eigenvalue E n (k) has an associated eigenvector with a definite 
spin orientation. 

The Pauli equation with spin-orbit coupling has the hamiltonian 

H = - X -A x + Vr(x) + X -g ■ (VVHaO A (-iV x )) . (2.33) 
After Zak transform the corresponding fiber hamiltonian becomes 

H(k) = i(-iV„ + k) 2 + V r (y) + \o- ■ (VV r (y) A (-»V V + k)) (2.34) 

with periodic boundary conditions. H of p.33p is bounded from below. But 
otherwise the band structure is similar to the periodic Dirac operator. Proposi- 
tion l2.1l and Corollary |2 . 251 hold as stated. In the proof one only has to use the 
appropriate time-reversal operator, which is Tip — o- 2 ifj* in the 03-eigenbasis. 

2.4 Gap condition and smoothness 

Let us consider one of the periodic hamiltonians, -ff pcr , with fiber decompo- 
sition H(k). H pcr is adiabatically perturbed to H e . Very crudely the corre- 
sponding unitary groups should be close. To make such a notion quantitative 
a gap condition must be imposed. We denote by a(H) the spectrum of the 
self-adjoint operator H . 
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Gap condition: We distinguish a family of m physically relevant energy 
bands {Ej(k) , n<j<n+m— 1} = o~a(k). This family satisfies the gap 
condition if 

dist{cr (fc) , <r(H(k)) \ CTo(fc)} > 9 > for all k S M* . (2.35) 

We repeat that the gap condition is not a spectral condition for H pcr . Let us 

setP° = ©^r 1 ^ 

Under the gap condition the projector P{k) depends smoothly, in many 

cases even (real) analytically, on k. RanP°(fc) is spanned by the basis 
{<y2j(fc)}j=«....,n+m-i- If the m relevant energy bands have no crossings amongst 
each other, then ipj is necessarily an eigenvector of H (k) satisfying H(k)<pj (k) = 
Ej{k)ipj{k). But if there are band crossings, it can be convenient not to insist 
on <Pj(k) being an eigenvector of H(k). Thus, while P°(k) is unique, the span- 
ning basis is not. In applications it is of importance to know whether there is at 
least some choice of ifj (k), j = n, . . . , n+m — 1, such that they have a smooth 
fc-dependence. Locally, this can be achieved. However, since M* has the topol- 
ogy of a torus, a global extension might be impossible. In fact this will gener- 
ically happen for the magnetic Bloch hamiltonian, see |DN80| INo81| |Lys85| 
for examples. Somewhat surprisingly, a reasonably general answer has been 
provided only recently [Pan06j . For the case of the Bloch hamiltonian, ana- 
lyticity has been proved before in cases d = 1, m arbitrary, and d arbitrary, 
m = 1, see Nenciu |Nen831 INen9T] and Helffer, Sjostrand |HS89j . They rely 
on analytical techniques. In |Pan06j topological methods are developed. 

Proposition 2.3. In case of the non-magnetic Bloch hamiltonian let either 
d < 3, to £ N or d > A, m = 1. Then there exists a collection of smooth maps 
M. d 3 k i— > ifij(k) € L 2 (M), j = n, . . . , n + m — 1, with the following properties 

(i) the family {yj(fc)}j=n,...,n+m-i is orthonormal and spans the range of 
P°{k). 

(ii) each map is equivariant in the sense that 

ifj(k) = T(X)ip(k + A) for all k e M d , A E A, (2.36) 

where r(A) is multiplication by e lX ' v . The same property holds for the non- 
magnetic periodic Dirac operator and Pauli operator with spin-orbit coupling. 

Remark. The proof uses the first Chern class of the vector bundle whose 
fiber at k is the span of the family {(pj(k)}j= nt ... t n+ m -i i-e. RanP°(fc). To 
establish continuity, and thus smoothness, this first Chern class has to vanish, 
a property, which does not hold for a magnetic Bloch hamiltonian except for 
some particular energy bands. 

If e is small, excitations across the energy gap are difficult to achieve. 
More precisely to P° one can associate a projection operator TI £ such that for 
arbitrary £, £' € N, r € K+, it holds 
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||(1 - iF) e -^ Et 77^|| < c Lt ,{r)e'U\\ (2.37) 

for < t < e~ l t with suitable constants c^/(t) independent of e. In other 
words that the subspaces II £ H and (1 — II £ )H almost decouple, i.e. decouple 
at any prescribed level of precision and over any polynomial length of the time 
span under consideration. For the specific case of the Bloch hamiltonian more 
quantitative details on the decoupling are provided in Section [3] 

If the gap condition is not satisfied, the dynamical properties are much 
more model dependent. Firstly the gap condition can be violated in various 
ways. In our context, since H(k) has discrete spectrum, the violation comes 
through band crossings. The behavior close to a band crossing has to be 
studied separately [LT051 [BT05 . In other models the energy band sits at the 
bottom of the continuous spectrum of H{k) without gap |Teu02] . Then an 
assertion like Equation (|2.37[) holds only under a suitable restriction to small 
I, £', usually 1,1' = I or perhaps I = 2, 1' = 1. 

The inequality (|2.37p makes no assertion about the dynamics inside the 
almost invariant subspace II e TL. While there is a general theory available 
[PST03a , we prefer to discuss the examples separately in the subsequent 
sections. 

3 Nonmagnetic Bloch hamiltonians: Peierls substitution 
and geometric phase corrections 

We discuss in more detail the effective dynamics for the Schrodinger equation 
with a periodic potential. For concreteness we fix the spatial dimension to be 
3. Under Zak transform the nonmagnetic Bloch hamiltonian becomes 

U z (^{-iV x - A(sx)) 2 + V r {x) + ^sx^Uz 1 = H% (3.1) 

with 

H% = \{ - iV v + k — A{ieVl)f + V r (y) + 0(feV£) . (3.2) 

Here V£ is differcntation with respect to k and satisfying the y-dependent 
boundary conditions (|2.6p . H e z is a self-adjoint operator on L^(R 3 , H 2 (M)), 
compare with (|2.8p . 

In (|3.2p we observe that the external potentials couple the fibers. To em- 
phasize this feature we think of (|3.2p as being obtained through Weyl quan- 
tization from the operator valued function 

H (k,r) = l(-iV y + k-A(r)j 2 + V r (y) + (f>(r) (3.3) 

as defined on (r, k) G M 6 and acting on TC[ with fixed domain H 2 (M) , see 
|PST03bj for details. In this form one is reminded of the Weyl quantization 
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of the classical hamiltonian function h c \(q,p) = \p 2 + V(q) which yields the 
semiclassical hamiltonian 

H sc = ^(-ieV x ) 2 + V(x) (3.4) 

acting in L 2 (K 3 ). The analysis of (|3.4[) yields that on the time-scale e~ 1 t the 
wave packet dynamics governed by H sc well approximates the flow generated 
by h c \. In contrast, the adiabatic analysis deals with operator valued symbols, 
as in (|3.3[) . and has as a goal to establish that the dynamics decouples into 
almost invariant subspaces and to determine the approximate dynamics within 
each such subspace. 

To be specific, let us then fix throughout one band index n and let us 
assume that the band energy E n is nondegenerate and satisfies the gap condi- 
tion. Therefore we know that E n : M* — > R is smooth and we can choose the 
family of Bloch functions tp n (k), with H(k)ip n (k) = E n (k)ip n (k), such that 
tp n depends smoothly on k. For each £ £ N = {0, 1, . . .} there exists then an 
orthogonal projection 7T| on 7i T such that 

\\m , ni}\\ < c e e e+1 (3.5) 

for some constants eg. Integrating in time one concludes that the subspaces 
IIfTL T are almost invariant in the sense that 

||(1 - ni)e~^ e ' tH -nm < II^IRl + \t\) Cl e t+1 E~ 1 ' (3.6) 

for any £, I' £ N. Note that the adiabatic time scale, order e~ l , can have any 
power law increase, at the expense of choosing the order of the projection 
-ZTJ sufficiently large. Only for times of order e 1 / 6 one observes transitions 
away from the almost invariant subspace. The zeroth order projection is at- 
tached to the n-th band under consideration, while the higher orders are 
successively smaller corrections to ilp. To construct TTq one considers the 
projection onto the n-th band, \ip n (k)) ((p n (k)\, as an operator valued function 
with values in B(H.{). From it we obtain the minimally substituted projection 
\f n {k— A(r)))(tp n (k— A(r))\. Its Weyl quantization is e-close to the orthogonal 
projection IIq. 

The second task is to determine the approximate time-evolution on 7T|7i T . 
Since the subspace changes with e, it is more convenient to unitarily map 
nfTi T to an e-independent reference Hilbert space, which in our case is simply 
L 2 (M*). The dynamics on L 2 (M*) is governed by an effective hamiltonian. It 
is written down most easily in terms of a hamiltonian function hf : M* xl 3 -> 
R. h\ is a smooth function. We also may regard it as defined on R 3 x R 3 and 
r , *-periodic in the first argument. h\ admits the power series 

t 

3=0 
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with e-independent functions hj. The effective quantum hamiltonian is ob- 
tained from h\ through Weyl quantization. The index £ regulates the time 
scale over which the approximation is valid and the size of the allowed error. 

In |PST03b] we provide an iterative algorithm to compute hj . In practice 
only ho and h\ can be obtained, at best ft 2 under simplifying assumptions. 
While this may look very restrictive, it turns out that already hi yields novel 
physical effects as compared to ho- Even higher order corrections seem to be 
less significant. 

To lowest order one obtains 

h (k,r)=E n (k-A(r)) + <f>(r), (3.8) 

which Weyl-quantizes to 

W £ [h ] = E n (k - A(ieV k )) + 0(teV fc ) (3.9) 

acting on L 2 (M*), where iV k is the operator of differentiation with periodic 
boundary conditions. (The twisted boundary conditions appearing in (|3.2j) are 
absorbed into the unitary map of n^H to L 2 (M*).) In solid state physics the 
Weyl quantization (|3.9|) is referred to as Peierls substitution. (|3.8p . (|3 . 9[) have a 
familiar form. The periodic potential merely changes the kinetic energy ifc 2 of 
a free particle to E n {k). The main distinctive feature is the periodicity of the 
kinetic energy. For example, in presence of a linear potential cf>, <p{x) = — E ■ x, 
an electron, initially at rest, will start to accelerate along E but then turns 
back because of periodicity in k. 

To first order the effective hamiltonian reads 

ht{k, r) = (V^(r) - VE n {k) A B{r)) ■ A n (k) - B(r) ■ M n (k) , (3.10) 

with the kinetic wave number k = k — A(r). The coefficients A n and M. n are 
the geometric phases. They carry information on the Bloch functions ip n (k). 
A n is the Berry connection given through 

A n {k) = i(f n (k) , ^kfn{k))n f (3-11) 
and M. n is the Rammal- Wilkinson phase given trough 

M n {k) - \i{V k y n {k), A{H{k) - E n {k))V k <p n {k)) Hf . (3.12) 

The Bloch functions ip n are only determined up to a smooth phase a(fc), 
i.e. instead of ip n {k) one might as well use e~ la ^ (p n (k) with smooth a : M* — > 
R. Clearly A4 n is independent of the gauge field a. On the other hand, A is 
gauge-dependent while its curl 

J7 n = VA^„ (3.13) 
is gauge independent. From time-reversal one concludes that 
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fl n (-k) = -fl n (k). (3.14) 

In particular, in dimension d = 2 for the first Chern number of the Bloch 
vector bundle one obtains 

/ dkQ n {k) = Q, (3.15) 
Jm* 

For the magnetic Bloch hamiltonian, (|3.14|) is violated in general, see Sec- 
tion [3J The integral in (|3.15p can take only integer values (in the appropriate 
units) and the first Chern number may be different from zero. Physically this 
leads to the quantization of the Hall current |PST03bl ISN99] . 

We still owe the reader precise a statement on the error in the approxi- 
mation. At the moment we work in the representation space at precision level 
I = 1. Let H e s be the Weyl quantization of ho + eh\, see p.8j) and (|3.10D . 
There is then a unitary U e : n{U T -> L 2 (M*) such that for all tp EH T 

\\( e - iH ^ -U s *e- iH ' t!t U e )nii;\\ < c||V||(l + |r|)£ 2 (3.16) 

with \t\ < e t and some constant c independent of ||?/>||, t, and e. 



4 Magnetic Bloch hamiltonians: the Hofstadter butterfly 



We turn to a magnetic Bloch hamiltonian in the form (|1.4|l , in dimension d = 2 
and with a transverse constant magnetic field Bq. We want to explain how 
the limits Bq — > oo and Bq — > can be understood with adiabatic methods. 
As a remark, it is worthwhile to recall that, when the physical constants are 
restored, the dimensionless parameter Bq is given by 

B a = B ° S (4 1) 

2i:hc/e ' 

where S is the area of the fundamental cell of _T and Bq the strength of 
the magnetic field, both expressed in their dimensional units. Thus -Bo corre- 
sponds physically to the magnetic flux per unit cell divided by hc/e, as the 
fundamental quantum of magnetic flux. This section is based essentially on 
[FP06] , which elaborates on previous related results |Bel86|. IHS89j . 

Adiabatic limits are always related to separation of time-scales. In the 
present case, one indeed expects that as Bq — * oo the cyclotron motion induced 
by Bq is faster than the motion induced by Vr , while in the limit Bq —> the 
microscopic variations of the wave function induced by Vp are expected to be 
faster than the cyclotron motion. 

Let us focus first on the Landau regime Bq — * oo. In order to make quan- 
titative the previous claim, one introduces the operators 

L\ = ^7= (pi + \Bqx 2 ) , 

[L x ,L 2 ]=i\, (4.2) 

L ^ = ^{P2-\BqX X ) , 
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and the complementary pair of operators 
G i = W (Pi ~ l B o a*) . 

[G!,G 2 ]=^1, (4.3) 

G 2 = i(p 2 + ±B *i) , 

where the relative sign is chosen such that [L^ Gj] = 0, for i, j = 1, 2. 

If Vp = 0, then Hmb describes a harmonic oscillator, with eigenfunctions 
localized on a scale |-Bo| _1//2 ; this corresponds to the cyclotron motion in 
classical mechanics. Since [Gj,Hmb] = 0, the operators G\ and G2 describe 
conserved quantities, which semiclassically correspond to the coordinates of 
the center of the cyclotron motion. 

If Vr 7^ 0, but the energy scale ||Vr|| is smaller than the cyclotron en- 
ergy ss £?o, then the operators Gi have a non-trivial but slow dynamics. By 
introducing the adiabatic parameter 77 = 1/Bq the hamiltonian reads 

Hmb = ^-{L\ + L\) + V r (G 2 - VVL2, -G x + ^1) • (4.4) 

In view of the commutator [G 1; G 2 ] = zr^l , one can regard r] Hmb as the 
77-Weyl quantization (in the sense of the mapping (q,p) > (Gi,G 2 )) of the 
operator-valued symbol 

ftMB(g,p) = 5 ( L ? + L i) + V Vr (p - -9 + 07 £ O • (4-5) 

For each fixed (q,p) £ K 2 , hyiB(q,p) is an operator acting in the Hilbert space 
Hi = L 2 (K) corresponding to the fast degrees of freedom. If ||Vr||s(w) < 00 1 
then /imb(?)P) has purely discrete spectrum, with eigenvalues 

A n ,„(<z,p) = („ + !) + rjVrip, -q) + G(r? 3/2 ), n G N, 

as 77 I 0. The index n G N labels the Landau levels. For 77 small enough, each 
eigenvalue band is separated from the rest of the spectrum by an uniform gap. 
Thus we can apply space-adiabatic perturbation theory to show that the band 
corresponds to an almost-invariant subspace 77 nj)7 L 2 (R 2 ). Let us focus on a 
specific n € N. One can prove that the dynamics inside Rani7 JM) L 2 (R 2 ) is de- 
scribed by an effective hamiltonian, which at the first order of approximation 
in r\ reads 

&? = (»+i)+»?VHGi,-G a ). (4.6) 

The first term in (|4.6p is a multiple of the identity, and as such does not 
contribute to the dynamics as far as the expectation values of observables 
are concerned. Leading-order dynamics is thus described by the second term, 
which does not depend on the Landau level n G N. Since Vr is a biperi- 
odic function and (Gi, G 2 ) a canonical pair, the second term is a Harper-like 
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operator. The spectrum of such operators exhibit a complex fractal behav- 
ior (Hofstadter butterfly) sensitively depending on the diophantine properties 
of a = l^r (notice that Vr{Gi,G2) depends on a through the commutator 
[Gi,G2] = iB^ l \). The Cantor structure of the spectrum was proven first 
in [BS82 for the case Vr(x\ 1 X2) = A cos a; 1 + cos x 2 (Harper model), for a 
dense set of the parameter values. Later Helffer and Sjostrand accomplished a 
detailed semiclassical analysis of the Harper operator H S89] . As a final step 
the Cantor spectrum has been proven by Puig (A ^ 0, a Diophantine) |Pu04j . 
and by Avila and Jitomirskaya [AJ05] for all the conjectured values of the 
parameters: A ^ 0, a irrational (the Ten Martini conjecture, as baptized by 
B. Simon). 

Secondly we turn to the opposite limit Bq —> 0, where the slow part of the 
dynamics is still described by the magnetic momentum operators Lj = V BaLj 
(j = 1,2), with commutator of order 0{Bq). However the decomposition given 
by (|4.2[) and (|4.3[) is no longer convenient. 

Since Aq is a linear function, Ao(ex) = t^eBq A x, the slow variation limit 
e — ► agrees with the weak field limit Bq — > 0. We then pose e = Bo and 
we regard Hmb in (|1-4|) as an adiabatic perturbation of the periodic hamil- 
tonian (12. 3|) . Thus we are reduced to the situation described in Section 3: to 
each isolated Bloch band of the unperturbed hamiltonian there corresponds a 
subspace 7T nie L 2 (R 2 ) which is approximately invariant under the dynamics as 
e I 0. The dynamics inside this subspace is described by Peierls substitution 
(|3.9|) . which now reads 

W £ [h ] = E n (k- ie 3 A(ieV fe )), (4.7) 

as an operator acting in L 2 (T 2 ,dk). Here Bq = (0,0, e) and e, is the unit 
vector in the i-th direction. 

Formula (|4.7p shows that the leading order effective hamiltonian depends 
only on the operators (K-l,K<2) = K. 

K = k-^e 3 A (ieV k ), 

which roughly speaking are the Fourier transform of the pair (L 1 ,L 2 ), and 
not on the complementary pair of operators. The same property holds true for 
the effective hamiltonian at any order of approximation £ £ N, see [FP06 , 
with important consequences on the splitting of magnetic subbands at small 
but finite Bq. 

An operator in the form (14. 7|) . shortly written E n (Ki, K2), is isospectral 
to an Harper-like operator, namely E n (G\ : G2) acting in i 2 (R). Indeed the 
first numerical evidence of the butterfly-like Cantor structure of the spectrum 
of Harper-like operators appeared when Hofstadter investigated the spectrum 
of cos K\ + cosK 2 as a function of e |Ho76j . On the other side, an operator 
of the form E n (Ki, K 2 ) is not unitarily equivalent to the Harper operator 
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E n (G\, G2). The important geometric and physical consequences of this fact 
are developed in [FP06j . 

Having explained the two extreme cases, -Bo — > and B a — > 00, the reader 
may wonder about the intermediate values of the magnetic field, Bo 1. As 
explained already in Section 12.31 it is convenient to introduce the magnetic 
translations 

T a = eT^ (a) exp(iB a-G), ae T , 

see (I2.17[) and (|2.19|) . If Bq satisfy Assumption 2, then {T a } is a commutative 
group, thus leading to the magnetic Zak transform (|2.20|) . Hmb is then a 
fibered operator over the magnetic Bloch momentum k G T 2 . At each k the 
spectrum of Hmb(k) is pure point and the corresponding eigenvalues are 
the magnetic Bloch bands. 

In view of this structure, one might argue that the adiabatic perturbation 
of the hamiltonian which includes, on top of the constant magnetic field Bq, 
a slowly varying magnetic potential A(ex) as in (|1.5[) can be treated with the 
methods of Section 3. There is however one crucial element missing. Indeed 
one can still associate to each magnetic Bloch band , isolated from the 
rest of the spectrum, an almost- invariant subspace Ran7T,f . On the other 
side the construction of the effective hamiltonian relies on smoothness which 
may be impeded for topological reasons. Indeed the analogue of Proposition 
12.31 is generically false for magnetic Bloch hamiltonians, as well-understood 
UDN80, N08I, Lys85 . In geometric terminology this fact is rephrased by saying 
that the magnetic Bloch bundle is generically non-trivial (in technical sense). 
This important fact has sometimes been overlooked. For example, Assump- 
tion B in UDGR04J is equivalent to the triviality of the magnetic Bloch bundle. 
Under this assumption the magnetic case is already covered by the results in 
|PST03b| . Thus the problem of adiabatic perturbation of a generic magnetic 
Bloch hamiltonian appears to be an open, in our view challenging, problem 
for the future. 



5 Piezoelectricity 

In the year 1880 the brothers Jacques and Pierre Curie discovered that some 
crystalline solids (like quartz, tourmaline, topaz, . . . ) exhibit a macroscopic 
polarization if the sample is strained. 

It turns out that also this effect can be understood in the framework 
of adiabatically perturbed periodic hamiltonians, cf. PSpT06, Lei05 . The 
perturbation is now slowly in time, 

H PE (t) = -^A x + V r(£t) (x,et). (5.1) 

If the potential Vp(x, et) has no center of inversion, i.e. there is no point with 
respect to which the potential has space-reflection symmetry, then the slow 
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variation of the periodic potential is expected to generate a non-zero current 
and can be shown to do so for particular examples A HI. !!)!)" . By translation 
invariance this current if averaged over a unit cell is everywhere the same and 
we denote the average current by J £ it) . For the following discussion we assume 
that Vr varies only for times in the finite interval [0, T] . Integrating the current 
per volume over the relevant time interval yields the average polarization, 



AP £ = [ dtJ £ (t). 
Jo 



In this section we discuss results that relate the current J £ {t) directly to the 
quantum mechanics of non-interacting particles governed by the hamiltonian 
(|5.ip , without the detour via the semiclassical model. For this we need to solve 
the Schrodinger equation with initial state p(0) — P(0) being the spectral pro- 
jection of -ffpE(O) below the Fermi energy E(0). Since the piezo effect occurs 
only for insulators, we can assume that E(0) lies in a gap of the spectrum 
of Hpe(0) and, in order to simplify the discussion, we also assume that this 
gap does not close in the course of time. Hence there is a continuous function 
E : [0, T] -> K such that E(t) lies in a gap of H PE (t) for all t. The state at 
time t is given by 

P £ {t) = u £ (t,o)p(o)u £ (t,oy, 

where the unitary propagator U £ (t, 0) is the solution of the time-dependent 
Schrodinger equation 

i£—U s (t,0)=H FE (t)U s (t,0) with [/ e (0,0) = l. (5.2) 
With the current operator given by 

f:=-[H(t),x]=--V x , (5.3) 

e e 

and the trace per volume defined as 

T(A) := lim J-R e Tr(l A „^) , (5.4) 

with 1a„ being the characteristic function of a 3-dimensional box A n with 
finite volume \A n \, the average current in the state p £ (t) is 

J°(t)=T(p%t)f). 

Finally the average polarization is 



AP £ = f dt T{p £ {t)J £ 
Jo 



(5.5) 



which is the main quantity of physical interest. The given framework allows 
us to describe the macroscopic polarization of a solid by a pure bulk property, 
i.e. independently of the shape of the sample. 
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In the simplest but most important case (see Paragraph (ii) in Section 1 
for a discussion of the model), the periodic potential Vr(x,st) is periodic 
with respect to a time- independent lattice r. For this case King-Smith and 
Vanderbilt KSV93 derived a formula for AP based on linear response theory, 
which turned out to make accurate predictions for the polarization of many 
materials. Their by now widely applied formula reads 

n c 

AP = —- J2 / dk (Mk, T) - A n (k, 0)) , (5.6) 

where the sum runs over all the occupied Bloch bands and A n (k,t) is the 
Berry connection coefficient for the n-th Bloch band at time t E M, 

An(k,t) = i((p n (k,t),V k (p n (k,t)) L 2( M} . 

Although A n depends on the choice of the Bloch function tp n , the average 
polarization (|5.6|) defines a gauge invariant quantity, i.e. it is independent of 
the choice of Bloch functions. 

In |PSpT06] we show that AP e defined in (|5.5p approaches AP as given 
by the King-Smith and Vanderbilt formula (|5.6p with errors smaller than any 
power of e, whenever the latter is well defined. More precisely we show that 
under suitable technical conditions on Vr{t) the average polarization is well 
defined and that for any N £ N 

£P° = -J— Fdt [ dkG(k,t)+0(e N ), (5.7) 
\ Z7T ) Jo Jm» 

where 

G(k,t) :=-itr (P(k,t)[d t P(k,t), V fc P(fc,t)]) , (5.8) 

and P(k, t) is the Bloch-Floquet fiber decomposition of the spectral projector 
P(t) = l(-oo.B(t)] (HpE(t))- Whenever all Bloch bands within RanP(fc,i) are 
isolated, the explicit term in (15. 7|) agrees with (|5.6|) . Note however that (|5.7j) 
is more general, since it can be applied also to situations where band crossings 
occur within the set of occupied bands. 

From the point of view of adiabatic approximation, this result is actually 
quite simple, since one just needs the standard time-adiabatic theory. At time 
t = the state p(0) is just the projection P(0) onto the subspace of the iso- 
lated group of occupied bands. Since these bands remain isolated during time 
evolution, this subspace is adiabatically preserved according to the original 
adiabatic theorem of Kato Ka50j, i.e. 

p z(t) = P(t) + 0(e), 

and one can compute the higher order corrections to p 6 (t) using the higher 
order time-adiabatic approximation due to Nenciu [Ne93j . In order to get 
explicit results, one has to do the adiabatic approximation for each fixed 
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k e M* separately. This is possible since Hps(t,k) is still fibered in k, due 
to translation invariance with respect to a time-independent lattice. However, 
since we need to differentiate with respect to k in order to compute the current, 
as suggested by formula (|5.8p . the expansion needs to be done uniformly on 
spaces of suitable equivariant functions. This makes the proof more technical 
than expected at first sight. 

Alternatively one can derive also for HpE(t) the semiclassical equations of 
motion including first order corrections: 

q = V k E n (k,t)-eO n (k,t), 



And again averaging the velocity over the first Brillouin zone yields the correct 
quantum mechanical average current that is the contribution from the n-th 
band. 

Note the striking similarity between the semiclassical corrections in (|5.8p 
and the electromagnetic field. If we define the geometric vector potential 

A n {k,t) = i(<p n (k,t), VktPn(k,t)) L z( M ), 

and the geometric scalar potential 

<f> n (k,t) = -i(ip n (k,t),d t (p„(k,t)} L 2 (M - ) , 

in terms of the Bloch function (p n (k, t) of some isolated band, then in complete 
analogy to the electromagnetic fields we have 

Q n {k,t)=-d t A n (k,t)-V k <f> n (k,t), (5.10) 

and 

fl»(M)=VjA^(M). (5-11) 

Time-dependent deformations of a crystal generically also lead to a time- 
dependent periodicity lattice r(t), see (15.11) . This more general situation is 
considered in [Lci05 , LP06J . Now the lattice momentum k is no longer a con- 
served quantity and the full space-adiabatic perturbation theory is required 
in order to compute the corresponding piezoelectric current. As a result an 
additional term appears in the semiclassical equations of motion, reflecting 
the deformation of the lattice of periodicity. 
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